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Abstract 

We outline the twist-two analysis of deeply virtual Compton scattering 
Q . (DVCS) within the conformal partial wave expansion of the amplitude, 

CN \ represented as a Mellin-Barnes integral. The complete next-to-leading 

order results, including evolution, are obtained in the MS and a con- 
Q . formal factorization scheme. Within the latter, exploiting conformal 

symmetry, the radiative corrections are evaluated up to next-to-next- 
^ ■ to-leading order. Using a new proposed parameterization for GPDs, 

we study the convergence of perturbation theory and demonstrate for 
r- 1 ! HI and ZEUS measurements that our formalism is suitable for a fitting 

Oh' procedure of DVCS observables. We comment on a recent claim of a 

O breakdown of collinear factorization and show that a Regge-inspired 

Q 2 scaling law is ruled out by small xbj DVCS data. 

j-C ■ 

\ 1 Introduction 

The proton structure has been widely explored in inclusive measurements, mainly in deeply in- 
t^J- \ elastic lepton-proton scattering (DIS). Here the scattering essentially occurs due to the exchange 

of a virtual boson (photon) between lepton and a single parton, and so one can access parton 
distribution functions (PDFs). These universal, however, convention-dependent functions q a {x) 
are interpreted as probabilities that partons of certain flavour a will be found with given longi- 
j^. ' tudinal momentum fraction x. Since the PDFs are naturally defined in a translation invariant 

manner, they do not carry information about the transversal distribution of partons. Some in- 
J> \ formation about transversal degrees of freedom can be obtained from elastic lepton-proton scat- 

\Ta ; tering. Namely, the electromagnetic form factors F\^{t) are Fourier transforms of the electric 

and magnetic charge distribution in nucleon, and can be, e.g., in the infinite momentum frame, 



interpreted as probability that partons are found at some transversal distance b from the center-of- 
mass. However, one should not assume that a realistic probability distribution of partons, given 
by a two-variable function q{x, b) is simply a direct product, i.e., q(x) ® q(h), of two probabil- 
ity functions. Rather it is anticipated that longitudinal and transversal degrees of freedom have 
a cross talk, e.g., as x gets bigger partons carry more of the nucleon longitudinal momentum 
and are expected to be closer to the proton center, and thus the b dependence in q(x, b) should 
become narrower with increasing x. 



' Based on talks presented by the authors at Elastic and Diffractive Scattering 2007, DESY Hamburg, Germany; 
Exclusive Reactions at High Momentum Transfer, JLab Newport News, US; The 6th Circum-Pan-Pacific Symposium 
on High Energy Spin Physics, UBC Vancouver, Canada; XII Workshop on High Energy Spin Physics, JINR Dubna, 
Russia; and New Trends in High-Energy PHYSICS, Yalta, Ukraine. 




Fig. 1: a) DVCS. b) Leading-order perturbative contribution to DVCS. c) Factorization (to all orders in a s ) on an 
example of special parity-even, helicity conserving contribution a Tt, Eq. d 1 Ob . 

The three dimensional distribution of partons in the nucleon can be addressed within more 
general objects, called generalized parton distributions (GPDs) [1-3]. Such distributions can 
be revealed by analyzing hard exclusive leptoproduction of mesons or photon. The theoretical 
description of the former processes is perhaps more problematic, but they offer a direct view 
into individual flavour GPDs. To the latter one the deeply virtual Compton scattering (DVCS) 
process contributes, where one photon has a large virtuality. DVCS is theoretically considered as 
the cleanest probe of GPDs, however, here only certain flavour combinations of GPDs appear. 

The non-forward Compton scattering process is described by the Compton tensor 



where q = (qi + q2)l% P = P\ + P2 and A = P% — P\. The generalized Bjorken limit 
corresponds to Q 2 = — q 2 — > 00 with the scaling variables 



and the momentum transfer squared A 2 being fixed. Note that in the forward case, i.e., A — > 0, 
the hadronic DIS tensor is related to the forward Compton scattering tensor by the optical 
theorem 



where p = P\ = P2 and £ — > XBj ■ 

In DVCS - Fig. [T^ - the virtuality of the incoming photon Q 2 = —qf is large while the 
final photon is on-shell. The skewness parameter 77 and the Bjorken-like scaling parameter £ 
are then equal to twist-two accuracy, i.e., 77 = £ + 0(1/ Q 2 ). In the generalized Bjorken limit, 
similarly as for DIS structure functions, the amplitude factorizes [4] into long- and short-distance 
contributions (FigfTJ) and c): short-distance physics is perturbatively calculable Compton scat- 
tering on a parton, while long-distance physics is encoded in a non-perturbative amplitude for a 
parton being emitted and later reabsorbed by the nucleon. The latter amplitude is called GPD. 
The factorization of long- and short-distance contributions is based on the reshuffling of divergent 
contributions from the perturbative calculable hard-scattering part to the universal, i.e., process 




(1) 



= ZmT^q, P = 2p, A = 0)/(2vr) 



(3) 
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independent GPDs. Although GPDs are universal they depend on the scheme convention and also 
on the order of the approximation to which one is able to perform the perturbative calculation. 

The GPDs with even parity, considered here, are defined as 



z+=0,z I =0 



G F(x,r,,A 2 ) = ±. f'^-j* p+z -(P 2 \Gi»(-z)G a +{z)\P l ) 

J Z7T z+=0, z x =0 

and similarly for odd parity. Furthermore, it is convenient to make a decomposition, 



(4) 
(5) 



ap = ujP^ujP,) aH u (P 2 ) ia +l/ u(Pi ) A„ aE 
P+ 2MP+ 

into helicity conserving and non-conserving generalized form factors. The Compton tensor (Q]) 
is analogously decomposed into Compton form factors with corresponding parity and helicity 
properties. In the forward limit (A — > 0) GPDs reduce to PDFs. Together with sum rules, e.g., 
relating GPDs to electromagnetic form factors, this provides constraints that are important for 
GPD modelling. But as they do not constrain much the skewness (77) dependence, modelling is 
still a difficult problem. One guiding feature is a polynomiality property of GPDs: n-th (Mellin) 
moment of GPD is even polynomial in r\ of order n or n ± 1. Similar polynomiality will be 
obeyed also by conformal moments which are just linear combination of Mellin moments. The 
usefulness of the GPDs has also been widely realized in connection with the spin problem, since 
they encode the angular momentum carried by the individual parton species, as explicated by the 
Ji's sum rule [3]. For recent detailed account of GPDs and their properties, we refer to [5]. 

In the following we first briefly comment on the work in Ref. [6], which claims that 
collinear factorization breaks down. We then proceed with the outline of the perturbative QCD 
approach to DVCS that is based on the conformal partial wave expansion, represented as Mellin- 
Barnes integral and described in detail in Refs. [7,8]. Here, as in Refs. [8-10], we concentrate on 
the dominant Compton form factor TC corresponding to parity-even helicity-conserving GPD H. 
For simplicity, we present only the results for the singlet part which is dominant for the kinemat- 
ics of collider experiments. We shortly discuss our phenomenological analysis of HI and ZEUS 
DVCS data. Finally, we conclude. 



2 Does collinear factorization break down? 

Here we would like to comment on a recent claim of an observed breakdown of collinear factor- 
ization [6], which the authors consider as a general feature of deeply exclusive leptoproduction. 
This claim originated from a divergent result that was found by convoluting a model GPD with 
the hard-scattering part in the handbag approximation, cf. Fig.QJ). In our opinion the GPD model 
employed by Ref. [6] suffers from ultraviolet and/or infrared divergencies. A matching proce- 
dure with the collinear factorization approach is impossible, since the model does not respect the 
ultraviolet behavior of QCD. 

• Divergences in a non-perturbative model that can not be cured within the factorization 
approach have to be removed by hand, e.g., by an effective cut off, before the model can 
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be utilized. Convoluting a model GPD, even a finite one, with the hard-scattering part, and 
obtaining undefined result does not in itself constitute a proof of a breakdown of collinear 
factorization. 

We recall that the collinear factorization approach has been worked out for deeply exclusive 
leptoproduction to NLO [11, 12], including evolution, and its quantitative features have been 
extensively studied. Below we report about our own investigations of DVCS beyond this order. 

We find it interesting that the analysis of Ref. [6] resulted in an infinite DVCS amplitude. 
Since it was claimed that this is a necessary consequence of Regge behavior, we would like to 
offer our opinion about the origin of these divergences. Let us remind that the field-theoretical 
definition (@| of GPDs states that they are generalized functions in the mathematical sense rather 
then regular ones [1]. To shorten the discussion and without loos of generality, we will employ 
this fact for the evaluation of the convolution integral [13]. 

The non-perturbative dynamics of the GPD utilized in Ref. [6] is described by a parton 
model. It is supposed that the collective spectator has leading Regge behavior, i.e., < a < 1, 
and interacts point-likely with the active spin-zero constituents. This reduction of the spin content 
induces a factor x which is taken outside of the GPD and leads to a non-standard representation 
[14]. This causes a mathematical paradox for "good" GPDs. Namely, the inverse, i.e., 1/x, 
moment of a GPD is defined, contrarily to what one expects on general grounds, while in the 
forward limit the inverse moment of the resulting PDF is divergent for < a < 1. Obviously, 
to be consistent, a non-standard GPD must contain some "bad" terms. The solutions are (i.) to 
reject such a GPD model, (ii.) to put the "bad" part in the central region — rj < x < rj, where it 
serves as a counterterm, and (iii.) to define the inverse moment for the PDF. Our choice in Eq. 
© below corresponds to the last possibility and the illustration of thus obtained "good" GPDs. 
Furthermore, we stress that the point-like coupling of the constituents ties the high-energy and 
the ultraviolet behavior. The treatment of divergencies as part of modelling was not specified in 
Ref. [6] and, consequently, the resulting GPDs contain ambiguous terms. 

• We consider the GPD model [6] in its presented form as mathematically ill-defined. 

An elegant method to overcome all these obstacles is to consider the GPD model [6] as an 
analytic function of the parameter a. This allows removal of ultraviolet (or infrared) divergencies, 
except at discrete values of a, e.g., a = and a = 1. Possible version of the GPD model part 
that originates from the leading Regge behavior reads in the region — rj < x < 1 as follows 



Up to terms which live only in the central region — rj < x < rj, and Regge non-leading terms, 
this expression coincides with the result given in Ref. [6]. Note that we consider (•••)" as 
generalized functions, which might be expressed using conventional "+" definitions and some 
additional subtraction terms. The forward limit yields the PDF H(x,rj = 0) = Nx~ a (l — x) 

'in the model one can choose to have ultraviolet, infrared, or both kinds of divergencies, however, one can not 
get rid of all of them. After momentum integration the GPDs will always suffer from the same illness, which is not 
primarily related to the high-energy behavior, rather, it is the generic feature of the underlying simplifications. 



„, \ N X „, . fx + 7]\ 

H(x,ri) = -- eix + rj j— M 

2a rj \ 1 + r] J 



—a 



(7) 
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and its inverse moment, i.e., fL\dx — x) = l/a(a — 1), is denned. Armed with a clear 

mathematical language [13], we calculate the convolution with the hard scattering part and get a 
finite result, however, with an 'unusual' negative overall sign: 



no 



dx ■ 



2x 



it 



-H(x,t) 



Ntt 



cot(^)+O(£)]+O(£ ). (8) 



• The reported divergences in Ref. [6] do not show up in the leading Regge behavior. They 
arise from ill-defined subtractions, entirely related to non-leading Regge behavior. 

Our inspection revealed that the ambiguous part in the GPD model [6], affected by ultraviolet 
and/or infrared divergences, possesses non-leading Regge behavior. It shows up only in the 
central GPD region and causes the observed divergences in the DVCS amplitude. 

The fact that the GPD representation is not unique was pointed out in Refs. [14, 15]. We 
recall that numerous QCD model evaluations confirm within the standard GPD definition a finite 
and continuous behavior at the cross-over points x = ±r/, see Refs. [5,16] and references therein. 
In contrast to the representation used in Ref. [6], Regge behavior of standard GPDs is continuous 
at the cross-over points, e.g., for a < 1: 



V 



1 — 77 



x + 77 

1 + Tj 



X + 7/ 
1 + 7/ 



1-a 



1] 



1-a 



spin-0 partons 



spin- 1/2 partons 



for x —* 77 , 77 < x . 



(9) 



If one chooses the standard GPD representation, the s-channel view in Ref. [6] confirms the 
implementation of Regge behavior in common GPD models [16] and our approach. To see 
this, one simply convolutes the spectral function (9) of Ref. [6] and the models (19) and (21) of 
Ref. [17] with respect to the spectator mass A. Doing so one finds the r.h.s. in (©. 

Finally, in Ref. [6] it was suggested to utilize so-called Regge-exchange amplitudes for 
phenomenology. These amplitudes are not universal and possess a typical s a ^ ~ (Q 2 /£) a( ^ 
Regge behavior, which yields a softening of Bjorken scaling. Let us recall that Regge phe- 
nomenology has been developed out of S'-matrix theory, i.e., for processes with asymptotic states, 
e.g., q\ = g| = 0, in lack of a dynamical understanding of strong interaction phenomena. Hence, 
one might surmise that the implementation of Regge-behavior in a parton model within virtual 
constituents and point-like interaction may predict unrealistic scaling for off-shell amplitudes. A 
last objection arises when the claimed scaling behavior [6] is confronted with DVCS data in the 
high-energy region (see Sect. 5 below): 

• The generic (Q 2 /£) a ^ behavior of Regge-exchange amplitudes is ruled out by experi- 
ment. 

We emphasize that this statement is based on experimental data, presently available. Our inves- 
tigation, given below, will show where perturbation theory reveals itself and in which kinematics 
the high-energy limit spoils the Bjorken limit. It turns out that this is rather a universal feature, 
which is tied to the ultraviolet behavior in QCD. 
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3 Deeply virtual Compton scattering 

Besides DVCS, the Bethe-Heitler (BH) brehmstrahlung process contributes to the measured hard 
photon leptoproduction off a proton. The BH amplitude, known in leading order of the QED fine 
structure constant, is expressed in terms of the known electromagnetic form-factors. Generally, 
there are two types of DVCS experiments: collider experiments and fixed target experiments. 
The former usually provide information in the phase space 10~ 4 < £ < 1CP 1 and lGeV 2 < 
Q 2 < 100 GeV 2 , and such are HI and ZEUS experiments at HERA. The main observables here 
are total and differential DVCS cross sections, however, also the measurement of beam charge 
asymmetry is feasible. For the fixed target experiments, such as Hall A, Hall B (CLAS) at JLAB, 
and HERMES at DESY, the DVCS-BH interference term can be more easily accessed via single 
beam, target spin and beam charge (HERMES) asymmetries, while the investigated phase space 
covers the so-called valence quark region, i.e., 0-05 < £ < 0.3 within 1 GeV 2 < Q 2 < 10 GeV 2 . 

One can express TC as a convolution (FigQ]:) over the longitudinal momentum fraction x 

a H(H, A 2 , Q 2 ) = J dx °C(x, e, Q 2 //x 2 ) a H(x, r, = £, A 2 , ^ 2 ) , (10) 

where fi 2 is a factorization scale that separates short- and long-distance dynamics and is often 
taken as \j? = Q 2 . Here the index a € {NS,S(S, G)} denotes either non-singlet or singlet parts, 
where to latter both quarks (S) and gluons (G) contribute. 

The coefficient functions C a are perturbative quantities which describe qj* — > qj and 
91* —> 91 subprocesses. The well known leading-order (LO) contribution to C a is actually a 
pure QED process (Fig. Q})). The next-to-leading order (NLO) contribution — the first order in 
a s — has been calculated by various groups [11]. Obviously, to stabilize the perturbation series 
and investigate its convergence one needs the second order in a s , i.e., next-to-next-to-leading 
order (NNLO) contributions. The importance of NNLO in singlet case is amplified by the fact 
that at LO photons scatter only off charged partons, whereas gluons start contributing at NLO. 

The GPDs a H are intrinsically non-perturbative quantities whose form at some initial scale 
Qo has to be deduced by some non-perturbative methods (lattice calculation, fit to data, etc.). The 
evolution to the factorization scale of interest is governed by perturbation theory, cf. [8, 10] for LO 
examples. The anomalous dimensions of non-diagonal operators were calculated up to NLO [18]. 
Still, evolution at NLO is numerically not easy to implement, and has been investigated beyond 
NLO only recently, using the procedure explained below [8]. 

Instead of using the convolution (TTOb one can equivalently use the sum over the conformal 
moments, and the singlet contribution then takes the form 

oo 

^(e, A 2 , Q 2 ) = 2 C'-'CjiQ 2 /^, a s (n)) = £, A 2 , ^ 2 ) , (11) 

3=0 

where Cj = (pCj, Cj) and Hj = (^Hj, G Hj) are conformal moments. They are analogous 
to common Mellin moments used in DIS but the integral kernel x- 7 is replaced by Gegenbauer 
polynomials C^ 2 {x) and C^ 2 {x), which are solutions of LO evolution equations for quarks 
and gluons, respectively. Unfortunately, the series (ITTb only converges in the unphysical region. 
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Hence, it is necessary to resum this series, e.g., by means of the Mellin-Barnes integral [8] 



a H{£,A 2 ,Q 2 ) = - 
1% 



i + tan ( — 



(12) 

The advantages of using conformal moments, i.e., Mellin-Barnes representation, are man- 
ifold: it enables a simple inclusion of evolution, it allows for an efficient and stable numerical 
treatment, and it opens a new approach to modelling of GPDs. Finally, by making use of confor- 
mal operator product expansion (OPE) and known NNLO DIS results, it enables the assessment 
of NNLO contributions, to which we now turn. 



4 Conformal approach to DVCS beyond NLO 

Neither Wilson coefficients nor anomalous dimensions are calculated in non-forward kinematics 
at NNLO (only so-called quark bubble insertions were partly evaluated [19]). To access the 
NNLO of non-forward Compton scattering, we use the conformal approach, making it possible 
to calculate relevant objects using only diagonal results of forward Compton scattering, i.e., DIS. 

DVCS belongs to a class of two-photon processes (DIS, DVCS, two-photon production of 
hadronic states . . . ) calculable by means of the OPE, T^(q, P, A) — > Cj Oj, for which the use 
of generalized Bjorken kinematics and conformal symmetry enables a unified description. While 
massless QCD is conformaily invariant at tree level, this invariance is broken at the loop level 
since renormalization introduces a mass scale, leading to the running of the coupling constant 
(J3 / 0). Assuming the existence of a non-trivial fixed point a*, i.e., (3(a*) = 0, the conformal 
OPE (COPE) prediction for Wilson coefficients in general kinematics reads [20] 

C 3 {a s ) - cM)iF\ (5 + 2i + 7i (a*))/2 e)W - (3) 

For 7] = equation (fT3T > reduces to the DIS Wilson coefficients Cj — > Cj and thus fixes the 
normalization cj. The choice -q = £ corresponds to DVCS in the conformal limit. The anomalous 
dimensions governing the evolution are diagonal and the same as in DIS. 

For a general factorization scheme, e.g., the MS scheme, the conformal symmetry breaking 
occurs also due to the renormalization of the composite operators and causes the appearance 
of non-diagonal anomalous dimensions = bjklj + T^JP ■ This induces a mixing of both 
operators, i.e., GPDs, and Wilson coefficients under evolution. In the kinematical forward limit 
(rj = 0) the diagonal evolution equations are again obtained, i.e., the DIS case corresponds to the 
COPE result. For DVCS, evaluated in the MS scheme, there appear also conformal symmetry 
breaking terms which are not proportional to f3, i.e., the non-diagonal terms survive. 

The non-diagonal terms of anomalous dimensions encountered in MS scheme can be re- 
moved by a finite renormalization [20], i.e, by a specific choice of the factorization scheme 

£yMs qMs = £fHs B B -i qMS _ ^CSqCS 



In this new scheme, called conformal subtraction (CS) scheme, all non-diagonal terms are "pushed" 



to the [3 proportional part 7^ s = Sjkjk+P/d^-jk- Furthermore, since there is an ambiguity in 
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Fig. 2: Relative NLO and NNLO corrections {L9} in the CS scheme (A 2 = 0.25GeV 2 , input scale Q 2 , = 2.5GeV 2 ). 
Thick [thin] lines denote "hard" ["soft"] gluon scenarios: Ng = 0.4, qg(0) = «e(0) + 0.05 [Ng = 0.3, qg(0) = 
cte(0) - 0.02]. 

MS — > CS rotation matrix, by judicious choice SB one can "push" mixing under evolution to 
NNLO. Hence, in CS scheme, which we are using, the unknown correction Ajk starts at NNLO 
and it can be additionally suppressed by the choice of an appropriate initial condition. Finally, 
we express our result in CS scheme as 



C 



cs 



J]C fc (a s (Q))Pexp 

k=j 



du! \ (3 
M L 9 



(15) 



with Ck(a s {Q)) obtained from the r] = £ limit of Eq. (fT3l ) and using Cj(a s ). As stated above 
the Akj mixing term appears at NNLO and is neglected. We take Cj and calculated to NNLO 
order from Refs. [21], and obtain the DVCS prediction to NNLO in the CS scheme. 



5 Results 

We have used the formalism described in the preceding sections to investigate the size of NNLO 
corrections to non-singlet [9] and singlet Compton form factors [10], to obtain complete (non- 
diagonal evolution included) MS NLO predictions [8], and to perform fits, in both MS and CS 
schemes, to DVCS and DIS data and to extract information about GPDs [8]. 

One can use a simple Regge-inspired ansatz for GPDs at the input scale 



ff;(f7,A 2 ,Q 2 o 



N^Fi;(A 2 )B(l+j-ax(0),8) 
N^F G (A 2 )B(l+j-a G (0),6) 



with (p a is a flavour dependent integer) 



a a (A 2 ) = a a (0) + 0.15A 2 , F a (A 2 



j + a a (0) 
j + a a (A 2 ) 



+ 0(7/2 



A 2 



Mf + jM^ 2 



(16) 



(17) 



Here a a are effective Regge-trajectories and F a at j = are partonic form factors, which are 
modelled as product of a monopole factor, arising from the effective Regge exchange, and an 
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impact form factor. The A 2 -dependence of higher moments could be easily pinned down from 
realistic lattice measurements, see Ref. [22] and references therein for the present status of sim- 
ulations in the flavour sector. The ansatz is based on modelling the t-channel contributions of 
the GPDs in terms of SO(3) partial waves [23] at a low input scale and then using analyticity, 
implicitly tied to Lorentz symmetry (polynomiality of GPD moments), to extend the GPD sup- 
port [7]. Our assumption is that at a rather low input scale, i.e., photon virtuality, Regge behavior 
is present in the GPD moments and can be effectively modelled by poles in the complex confor- 
mal spin plane. Doing that in such a way yields GPDs that behave continuously at the cross-over 
point and can so be smoothly incorporated in the factorization approach. We emphasize that this 
is compatible with the s-channel view [6], if it would be presented in an appropriate representa- 
tion, cf. Eq. (O. In the forward case (A = 0) ansatz (fTBT ) is equivalent to the standard building 
blocks for PDFs: 

E(x) = N^x- a *W (1 -x) 7 , G(x) = N' G x- aG{0) {l-xf. (18) 

We relied in our ansatz ([T6l ) only on the leading SO(3) partial wave, which can be for rj < 0.3 
safely approximated by a constant. The work on a more flexible r/-dependent ansatz, i.e., the 
model dependent resummation of SO(3) partial waves is in progress. 

We have performed the analysis of radiative corrections with generic parameters and made 
fits of relevant parameters N%, as(0), Mq, Nq, ac(0), Mq (since of j ~ dominance 
can be safely set to zero). We introduce now the quantities that we utilize as a measures of the 
scheme dependence and, foremostly, as indicators for the convergence of the perturbation series. 
It is natural to employ for this purpose the ratios of Compton form factors, i.e., the corresponding 
modulus and phase difference, at order N p LO to those at order N p_1 LO, where P = {0, 1, 2} 
stands for LO, NLO, and NNLO order, respectively: 



5 K 



yN p LO \ 

<^ = argl— ^— . (19) 



The phase differences are small, and we will not comment on them here further. The NLO cor- 
rections to the moduli in MS and CS schemes have a similar £-shape, where MS corrections are 
generally larger. The relative NLO and NNLO corrections in CS scheme are depicted in Fig. [2] 
From the left panel, showing corrections at the input scale, we realize that the large negative NLO 
corrections to the modulus in the 'hard' gluon scenario (thick dashed) are shrunk at NNLO to 
less than 10% (thick solid), in particular in the small £ region. In the 'soft' gluon case the NNLO 
corrections (thin solid) are ±5%. For £ ~ 0.5, the corrections are reduced only unessentially and 
are around 5% and 10% at both NLO and NNLO level. If evolution is switched on (right panel), 
our findings drastically change. For 5 ■ 10~ 2 < £ NNLO corrections are stabilized on the level of 
about 3% at Q 2 = 100 GeV 2 . But they start to grow with decreasing £ and reach at £ sa 10 -5 the 
20% level. As in DIS, this breakdown of perturbation theory at small £ in DVCS obviously stems 
from evolution and is thus universal, i.e., process independent. The large change of the scaling 
prediction within the considered order does not influence the quality of fits, and, in particular, the 
possibility of relating DVCS and DIS data. Hence, the problem of treatment or resummation of 
these large corrections is relevant primarily to our partonic interpretation of the nucleon content. 
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Fig. 3: Simultaneous fit to the DVCS and DIS data in the CS scheme to NNLO. Upper left panel DVCS cross section 
for Q 2 = 4 GeV 2 and W = 71 GeV (circles, dashed) as well as Q 2 = 8 GeV 2 and W = 82 GeV (squares, solid) 
[27] . Upper right panel DVCS cross section (| A 2 1 < 1 GeV 2 ) versus Q 2 for W = 82 GeV (HI, circles, dashed) and 
W = 89 GeV (ZEUS, triangles, dash-dotted) [28]. Lower left panel DVCS cross section versus W Q 2 = 4 GeV 2 
(HI, circles, dashed), Q 2 = 8 GeV 2 (HI, squares, solid), and Q 2 = 9.6 GeV 2 (ZEUS, triangles, dash-dotted). Lower 
right panel shows F 2 (x B j, Q 2 ) versus Q 2 forxsj = {8 ■ 1CT 3 ,3.2 ■ 1(T 3 , 1.3 ■ 1(T 3 , 5 ■ 1(T 4 } [29]. 



As long as we precisely define the treatment of the evolution operator, perturbative QCD can be 
employed as a tool for analyzing data also in the small £ region. However, it might be expected 
that in meson leptoproduction dominant perturbative corrections [12, 24] at small £ in the hard 
scattering amplitude are not connected only to the evolution and should be resummed [25]. 

The Mellin-Barnes integral approach offers the possibility for a fast and numerically stable 
analysis. Our numerical routine is designed for the purpose of fitting DVCS (and DIS) observ- 
ables and testing various GPD ansaetze. A fit example for measurements of the HI and ZEUS 
collaborations [26-29] is presented in Fig. [3] at NNLO. In these collider experiments the DVCS 
cross section can be accessed and is given, up to £ and (— A 2 )/4M 2 ~ 0.05 proportional correc- 
tions, at leading twist as 

^ WA ,, S ^^ |H|2 ( 5 ._^_, A2 , Q ,). 

For fixed £ we have the canonically expected scaling behavior 1/Q 4 , which is strongly modified 
by the resummation of log Q 2 terms. 

As one can realize from Fig. [3] the quality of a simultaneous DVCS and DIS fit within our 
ansatz (fT6l) is quite good, namely, x 2 /d.o.f. = 0.77. We only observed bad fits to LO accuracy, 
which originate from the fact that the normalization of the DVCS amplitude for fixed A 2 is not 
adjustable, since we took only the leading SO(3) partial wave. We already demonstrated that 
taking the next-to-leading SO(3) partial wave leads to good LO fits [30]. Note that this is nothing 



10 



else but a minimal version of the so-called dual GPD parameterization [23], used in Ref. [31]. Let 
us emphasize that our fits support the following functional form of CFFs in the small £ -region: 



where N(£, Q 2 ,t) and a eS (t, Q 2 ) effectively describe the change of the normalization and slope 
due to the predicted scaling violations, i.e., mainly arising from evolution. 

In contrast, the so-called Regge-exchange amplitude [6] yields a behavior which is rather 
similar to the Regge behavior of on-shell amplitudes (W 2 = s) 



Here the normalization factor iV'(£,A 2 ) can for internal consistency of the model [6] only 
weekly depend on £ ~ Q 2 /2W 2 and so it can be safely neglected in the small £ kinematics. 
This prediction has so far been compared only with exclusive leptoproduction data [32-34] in 
the valence quark region [6], which is not conclusive H. As spelled out above, the DVCS ampli- 
tude can be described to LO accuracy, where gluons are absent in the hard scattering part. Hence, 
for the sake of illustration we consider it legitimate to confront the Regge-exchange amplitude [6] 
within a pomeron Regge trajectory a(A 2 ) = 1 + 0.25A 2 < 1 with high-energy DVCS data. As 
one can immediately realize from Eqs. (1201 ) and (1221 for fixed W the cross section should not 
scale with Q 2 , which is in conflict with the observation, describable by a (Q 2 ) 1 - 54=F ° 09=F ° 04 
fit [27], as shown in the upper right panel in Fig. ©. Confronting Eq. (|22l with the few ex- 
perimental data points for fixed £ [27] also disfavors the claimed Regge scaling. We note that a 
BFKL evaluation of the DVCS amplitude indicates a rather intricate £ and Q 2 / A 2 dependencies 
in the high-energy limit [35], which seems not to support the conjectured Regge behavior [6], 
borrowed from on-shell amplitudes. 

Turning now to presentation of the GPDs resulting from the fits, we recall that Fourier 
transform of GPDs for rj = 0, 



can be interpreted in the infinite momentum frame as probability density [36], see Fig. |4Jd. The 
average transversal partem distance squared (b 2 ) is given by the GPD slope B = (b 2 )/4, shown 
in Fig. |4^. Although gluons are perturbatively suppressed in DVCS, one has a handle on them 
via the Q 2 evolution. The found value for (b 2 ) is compatible with the analysis for deeply exclu- 
sive photo and leptoproduction of J/ip [37]. The results confirm the picture, mentioned in the 

2 Note that from present fixed target data it is rather difficult to judge on any predicted Q 2 scaling behavior, since 
kinematical variables are strongly correlated for a large Q 2 lever arm and, moreover, to access the twist-two sector in 
exclusive meson leptoproduction one has to extract the cross section that arises from the exchange of a longitudinal 
polarized photon. Within small Q 2 lever arm and small t, DVCS [33] indicates a canonical scaling which might be 
logarithmical modified [8], and is also still consistent with Regge scaling within a rather low a(<C t ^) = 0.15 
parameter [5]; as in DIS, exact canonical scaling behavior perhaps takes place only in a very limited £ -region [8]. 




(21) 




(22) 




(23) 
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a) b) 
Fig. 4: a) Resulting GPD slope B = (b 2 )/4 at the input scale Q 2 = 4 GeV 2 and b) 3D picture of gluon GPD <23l . 



introduction, about the correlation of transversal and longitudinal degrees of freedom: partons 
with smaller momentum fractions are more decentralized. Intuitively, this small x partons might 
be somehow associated with the meson cloud that surrounds the proton center. 

6 Summary 

GPDs encode a unified description of the proton suucture and they are experimentally accessible 
via the hard leptoproduction of photon or mesons. We have shown that the representation of 
Compton form factors as Mellin-Barnes integrals offers a useful tool in analyzing DVCS: the 
inclusion of evolution is simple, numerical treatment is stable and fast. Although the motivation 
for this representation originated from manifest conformal symmetry at LO, we have shown here 
that this Mellin-Barnes integral representation can be used within the standard MS scheme be- 
yond LO. Such a representation can be straightforwardly obtained from the momentum fraction 
representation and, therefore, also other GPD related processes, e.g., the hard electroproduction 
of mesons, can be given in terms of Mellin-Barnes integrals. This opens a new road for the 
'global' analysis of experimental data within the perturbative GPD formalism to NLO accuracy. 
Furthermore, the use of conformal symmetry enables elegant approach to higher-order radia- 
tive corrections to the DVCS amplitude. We have shown that although NLO corrections can 
be sizable, and are strongly dependent on the gluonic input, the NNLO corrections are small to 
moderate, supporting perturbative framework of DVCS. The observed change in the scale depen- 
dence is not so conclusive: similarly as in DIS we encounter large NNLO effects for £ < 10~ 3 , 
which signal a breakdown of naive perturbation series expansion of the evolution operator. Nev- 
ertheless, this breakdown is universal and if we precisely define the treatment of this operator, 
perturbative QCD can be employed as a tool for analyzing data even in the small £ region. We 
expect that within increasing accuracy in the perturbative approximation the GPD formalism 
improves. 

Finally, fits to available DVCS and DIS data in collider kinematics work well in the 
collinear factorization approach beyond LO and give access to transversal distribution of par- 
tons. The lesson learned from the failure of LO fits can be extended to (oversimplified) GPD 
models used for phenomenology in fixed target kinematics. They do not possess a flexible skew- 
ness dependence parameterization and so the normalization of the amplitude for fixed t is mostly 
determined. Ironically, this non-flexibility originates from the implementation of PDF and form 
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factor constraints in the most convenient manner. This can be repaired in the Mellin-Barnes 
representation by a model dependent resummation of SO(3) partial waves or in the momentum 
fraction representation by a more flexible double distribution ansatz, in which also the correlation 
of t- and skewness dependence is improved [17]. 

Additionally, we have observed that the off-shell Regge-exchange amplitudes [6], having 
the generic s-dependence of high-energy on-shell amplitudes, are ruled out. Repairing this failure 
yields the well-known dilemma of Regge-phenomenology. Namely, when extended to off-shell 
processes it fails to describe high energy data and must be improved, thereby loosing its predic- 
tive power [38]. This loss of predictive power appears also in description of deeply exclusive 
leptoproduction data in the valence quark region within Regge inspired models, often outside of 
the accepted validity region {25 GeV 2 < W 2 , —t < lGeV 2 } in classical Regge phenomenol- 
ogy, see e.g. [39]. Certainly, these approaches offer an easy and economical possibility to fit data 
and to interpret them in terms of mesonic exchanges and numerous couplings (unknown form 
factors). In our opinion the GPD formalism offers a new perspective for our QCD understanding, 
including a handle on the dynamical sources of Regge behavior. In this formalism one extracts 
information about non-perturbative quantities from experimental observables and interprets it in 
terms of partonic degrees of freedom. We consider Regge-inspired interpretations of data [40] 
as valuable tool for building realistic GPD models. This gives one handle more on the challeng- 
ing task to understand hadron physics from QCD dynamics. Non-perturbative GPD aspects, in 
particular Regge behavior, might be studied within (partial) resummation of i-channel ladders 
(in the language of light-cone wave functions of higher partonic Fock states), e.g., Ref. [35,41], 
lattice simulations, and model building. This might lead to clues to improved understanding of 
hadron physics in terms of the underlying theory. 
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